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1. Introduction 

Ken Wilson introduced the exact renormalization group (ERG) to define continuum 
limits in quantum field theory, p] Using a theory with a finite ultraviolet (UV) cutoff, 
we can construct a continuum limit with the expense of keeping an infinite number 
of terms in the action. The idea was subsequently applied to perturbation theory by 
Joe Polchinski.[2j The purpose of this paper is to apply Polchinski's version of ERG to 
construct QED perturbatively. 

Before writing anything about gauge theories, we would like to clarify our 
terminology. The "action" we have mentioned above is a functional of field variables. 
The momenta of the fields are restricted to the range below a UV cutoff A. The 
interaction vertices of the action are obtained by integrating out the fluctuations of 
fields with momenta larger than A. Hence, the correlation functions calculated with 
this action contain the physics of all momentum scales. The action with a UV cutoff 
is often called a "Wilson action" in the literature, and we will follow this practice here 
even though an "action with a UV cutoff" is perhaps more descriptive. The reason is 
that we wish to avoid confusion with the "effective action" with an infrared (IR) cutoff 
A that consists of the 1PI parts of the vertices of the Wilson action. In the limit A — > 0, 
the "effective action" with an IR cutoff reduces to the standard effective action. ERG 
determines the UV cutoff dependence of the Wilson action or the IR cutoff dependence 
of the "effective action" . Except for this introduction, we will consider only the Wilson 
action in this paper. 

The construction of gauge theories in the framework of ERG has been discussed 
by various authors. To the author's knowledge the first person who applied ERG 
systematically to construct (non-abelian) gauge theories is Becchi. [3] (There is an even 
earlier work by Warr [I] , who treated ERG mainly as a method of regularization.) Becchi 
showed how to construct a BRST invariant Wilson action by introducing sources that 
generate BRST transformations. Alternatively gauge theories can be formulated for the 
effective action with an IR cutoff. This was initiated by Ellwanger [5] under the name 
of the "modified" Ward (or Slavnov- Taylor) identities. 

Whether we study a Wilson action with a UV cutoff or an effective action with an 
IR cutoff, the form of the Ward identities depends on the cutoff A explicitly. (This is 
the reason for the adjective "modified.") But ERG guarantees the A independence of 
the identities: once the identities are satisfied for some A, it is satisfied for any A. It is 
this property which simplifies the study of gauge symmetry within the ERG framework. 

The first explicit construction of QED using ERG was done in [6]. They computed 
the effective action with A = by solving ERG differential equations starting from a 
set of initial conditions at a very large cutoff A . They showed how to fine tune the 
initial conditions so that the effective action with A = satisfies the Ward identities. 
There, ERG is merely a device for loop calculations, and they did not take advantage of 
the simplicity of the Ward identities at large A. Except for the original work of Becchi, 
most later works seem to share the same short-coming: the Ward identities are always 
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checked at A = 0. 

The main purpose of this paper is to construct the Wilson action of QED by 
satisfying the Ward identities at large UV cutoff A. As advocated by Becchi [3], the 
Ward identities get simplified at large A, and to take advantage of this we need to 
parameterize the Wilson action in terms of its asymptotic behavior at large A. Such 
parameterization scheme has been introduced in refs. [TJ [8], and we apply it here to 



We organize this paper as follows. In sect. 2, we briefly review Polchinski's ERG 
differential equation by applying it to QED. In sect. 3 we parameterize the Wilson action 
in terms of its behaviors at large cutoff. This is an application of the scheme explained 
in [8]. (The reader unfamiliar with perturbative applications of ERG may find a casual 
reading of [8] helpful.) In sects. 4 and 5 we give a new derivation of the Ward identities 
for the Wilson action. The result we obtain is equivalent to what has been derived 
before by Becchi [3] and Ellwanger [5] and others for YM theories. Our starting point 
is the Ward identities in the continuum limit of QED. A crucial observation, which has 
not been emphasized or clearly understood in the previous literature on ERG, is that 
the continuum limit can be constructed using a Wilson action with any finite UV cutoff. 
(This is derived for the real scalar theory in appendix A.) Hence, we can rewrite the 
Ward identities in the continuum limit as the identities for the Wilson action. We will 
end up with an operator equation (1271) . an equality between two composite operators. 
In sect. 6, we explain how we can satisfy fTSTj) by fine tuning the parameters of the 
theory order by order in loop expansions. We give 1-loop calculations in sect. 7. After 
giving brief remarks about how to rewrite the Ward identity ([271) as BRST invariance 
in sect. 8, we conclude the paper in sect. 9. Some detailed calculations are given in 
appendix B. 

Before closing we mention some relevant past works. Besides [6], construction of 
QED using ERG has been done by others, and we mention only a couple here. In [9| [TO] 
the background field method was applied to scalar QED. More recently, a manifestly 
gauge invariant formulation has been constructed not only for QED [H] but also for 
general YM theories. [T2] 

Throughout the paper we use the euclidean metric. 

2. Polchinski's equation 

The Wilson action of QED is given as the sum of free and interaction parts: 



QED. 



S(A) = S free (A) + SUA) 



(1) 



where 



>Sfree(A) 



2J k 



k) Y{k/k) ( x i) kflK ) 




(2) 
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We use a sign convention so that the weight of functional integration is given by 
exp[5(A)] instead of the more usual exp[— S'(A)]. The momentum cutoff A is introduced 
via a positive cutoff function K(x), which is 1 for x 2 < 1 and decays rapidly (e.g., 
exponentially) for x 2 > 1. The free propagators are given by 

(AMM-k))^ = (v - (i - 0^) (3) 

mt(-p)) s (A) (4) 

Throughout the paper we consider only the connected part of the correlation functions. 
For notational simplicity we omit the overall factor of the delta function for momentum 
conservation. 

We can construct the continuum limit using a Wilson action 5(A) with a finite 
A. But this comes with a price: we must keep an increasing number of terms in the 
Wilson action as we go to higher orders in perturbation theory. The A dependence of 
the interaction part Si nt (A) is given by Polchinski's ERG differential equation [2]: 

8S int (A) affint(A) 

SA^k) 6A v (-k) ' 8A tl {k)8A v {-k) j 
1 





A(p/A)Sp 



im 



x { Tll—A^) ■ 5 int (A)-^- + - T ^_5 int (A)- 



5ip(—p) 8tp(p) 5ifj(—p) 8ijj(p) 



(5) 



A (t) = A ~ir K (t) (6) 



where 

A ( , _ . 

,AJ OA VA, 

is a non- negative function which is zero for p 2 < A 2 , has an appreciative value only for 
p 2 ~ A 2 , and decays rapidly as p 2 grows. The minus sign in front of the second integral 
of (JHJ) is due to the Fermi statistics. 

It is common to think that S(A) describes only the physics of low momentum p < A 
correctly. On a close inspection, however, we find that S(A) contains the physics of all 
momentum scales. To be precise, the renormalized two-point functions are given by 

{M-k)Mk)) = 1 " K fJ Ayl (V - (i - 0^) 

1 

+ K(k/A) 2 ^ A ^~ k ^ Av ^s(A) ( 7 ) 

{m ^- p)) = izJ^ + _i_(^ ( - ?) ) s(A) (g) 

where (• • -) s is the correlation function calculated with S. The renormalized higher 
point functions are more simply given by 

(Aui(^i) ' ■■A t t M (k M )ip(p 1 ) ■ ■ -ipiPN^i-qi) ■ ■ "0(-?7v)) 
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M i N y 

= 2 K^/Xj n K(pi/A)K(qi/A) 

• (A w (Ai) • • • A MM (k M )iP( Pl ) ■ ■ ■ 1>(pN)$(-qi) ■ ■ ■ $(-QN)) m ( 9 ) 

To make sense of the division by K(p/A), we assume that K(p/A) is strictly positive 
and that it vanishes only in the limit p 2 — > oo. As long as the Wilson action 5(A) 
satisfies (j5J), the right-hand sides of (J7f8|9l) are independent of the cutoff A. Since the 
results fl7|8|9|) may not be widely known, we derive analogous results for the real scalar 
field theory in appendix A. 

3. Parameterization 

Let us next discuss how to parameterize the solutions of (jSJ). We are only interested in 
those solutions corresponding to renormalized theories. (See 0, E] for a similar analysis 
for the 4 theory.) We evaluate the functional derivatives of S(A) at vanishing fields. 
The asymptotic behaviors at large A are obtained by expansions in powers of the electron 
mass and external momenta: 



5 2 S 



int. 



5^\a 2 (lnA/^A 2 



5A^k)5A v {-k) 

+ b 2 (lnA//i)m 2 + c 2 (hi A///) k 2 \ +d 2 (In A///) k± v H (10a) 



5ip(—p) Sip(p) 

t 5S *S 



Stp(-p- k) SA^k) 5ip(p) 
S 4 S 



a f (In A///) fi + bf (In A///) im-\ (106) 

= a 3 (hiA/ y u) 7A( + --- (10c) 



6A a (k 1 )SA/ 3 (k2)SA r (k 3 )5A 6 (k 4 ) 

= a 4 (In A///) (SapS-ys + 5 ay 5[3s + S a sSp 7 ) H (lOd) 

where an arbitrary momentum scale \x > is introduced to make — dimensionless. The 
vertical line with subscript is a reminder that the derivative is evaluated at vanishing 
fields. (The extraction of the asymptotic parts is similar to the T operation in [3].) 

For the renormalizable theories, the dotted parts are suppressed by negative powers 
of A, hence vanishing in the limit A//i — >• oo. The higher order derivatives have negative 
dimensions, and vanish also in the limit Aj\i — > oo. 

In solving ((SJ), the seven functions 

fo 2 (lnA///),c 2 (lnA/^),d 2 (lnA//i), 
ay (In A//x), 6/(ln A//x), a 3 (lnA//i), a 4 (ln A/fi) 

are ambiguous by additive constants. (Note that the quadratically divergent A 2 a 2 has 
no such ambiguity.) Therefore, in order to specify a unique solution of (jSJ), we must 
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introduce seven conditions. To begin, we can adopt the following three normalization 
conditions: 

qa(O) = o/(0) = 6,(0) = (12) 

(02(0) normalizes the gauge field, a,(0) the spinor fields, and 6/(0) the mass parameter 
m.) It is one of the main purposes of this paper to show how the Ward identities 
determine the remaining four constants 

6 2 (0),rf 2 (0),a 3 (0),a 4 (0) (13) 

in terms of the elementary charge e, which itself is introduced through the Ward 
identities. 



4. The Ward identities 

Let us recall the renormalized Ward identities in QED. They have two parts. The first 
part is the identity for the two-point function of the gauge field: 

!*„ (A,{-k)A v {k)) = ^ (14) 

The second part is the Ward identities for the higher-point functions: 

jkp (A^-fyA^h) ■ ■ ■ A^ M {k M )^{px) ■ ■■^(PN)^(-qi) ■ ■ •$(-qN) 

e N r / vi 

= p E [<An (*i) " " " 1>(Pi -*)••■)- (4* (h) ■ ■ -ip{—qi - k) ■ ■ ■)] (15) 

where e is interpreted as the elementary charge defined at scale \x. The Ward identities 
can be also written as the conservation law: 



(j^-tyA^ki) ■ ■■A flM (k M )'il)(pi) ■ ■■ip(p N )^(-qi) ■ • "0(-?jv) 

N 

= e £ [<A W (h) ■ ■ ■ i>{ V i ~ k) ■ ■ ■) - (A„ (h) ■ ■ ■ <K-g, - k) ■ ■ ■)} (16) 

i=i 

where is the charge current. The goal of this section is to find a concrete expression 
of J M in terms of the interaction action Si nt (A). 

We first examine (1T41) . Using ([7]), we can rewrite (1T4"|) as 

~k„ (A,(-k)A u (k)) s{A) = K(k/A)^ (17) 



Since 



(Ap(—k)A v (k)) s ^ 



K{h/A) f s»«- (i-e)^l ■ *«» + ( JhsMk)) i n> 



k~ V k 2 J \ \ . 1,» ( k ! / \ , 

(see appendix A for an analogous result for the real scalar theory), we can rewrite ([I 

as 
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Next, we examine (1151) . Using and 



5(A) 



/ ^int 

-^ 1 l/c 1 j---^lPiJ---^l-?iJ--- / 

5(A) 



££^4* • ■ • V»(pO ■ • ■ iK-fc) • ■ • ) (20) 



( again see appendix A), we can rewrite (fT5l as 



MV'V / 5(A) 

1 

X 



N 

£ 

n=l 



JT(tyA) njLi K( Pi /A)K( qi /A) 
(A^(ki) ■ ■ -ipipn - k) ■ ■ ■) - (A w (A;i)---0(-g n -A;)---) 



(21) 



To summarize so far, we can replace (I14|I15I) by (119)2ip . Note that (jT9l) is a special 
case M = 1. iV = of (Gill. Hence, (JSTJ) for arbitrary M, A is equivalent to fll4ll5D . 
We now define a vector current by 

This is a composite operator since its correlation functions 
(j M (-A:)A M1 (A; 1 )...V(p 1 )--^(-gi)---) 

= (^(-fc)^^) • • -V(pi) • ■ ■ <K-gi) ■ ■ -) 5(A) 

1 

nfl! K(h/A) nUK( Pj /A)K( qj /A) 
are independent of A. (See Appendix A for an explanation on composite operators.) To 
show this, we differentiate (jSJ) with respect to A^k) to obtain 
d_ 

dA' 

where T> is a linear differential operator defined by 



X i-rM TTfT, / a \ ttAT F77I /»u//_ /an ( 2 ^) 



- A— J M (-fc) = V ■ J^-k) (24) 



+ r- 



VM M (-/)^(/) 2M,(0M„(-I) 
A(p/A)Sp— ^- f-t—Ste'F- <] 



p 1 + im\5ip(—p) m d~i>(p) 

5ip(—p) m 6ip(p) 5ip(—p) d~vjj(p) J 
acting on any functional JF. (|24|) guarantees that the correlation functions (12"51) are 
independent of A. Thus, from (!2lT[23|) . we can rewrite the Ward identities as (|T6|) where 
is defined by (12"2"1) . (For a discussion of composite operators in 4 theory, see [8].) 
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5. Ward identity as an operator identity 

In this section we wish to rewrite the Ward identities ({TBI as a single operator equation. 
For this purpose we construct a composite operator $(—/&) whose correlation functions 
give the right-hand side of ( fT6l ): 

[<&{-k)A^ (ki) ■ ■ ■ A^ M (k M )ip(pi) ■ ■ • ^ipN)ip(-qi) ■ ■ ■ ii{-q N ) ] 

N 



e E (<4* ■ • • ^fe -*)••■>- ( fc i) • ■ ■ - fe) ■ ■ 



i=l 



(26) 



Once we construct $(— fc), the Ward identities (fTBl) are equivalent to the single operator 
equation: 

kM~k) = $(-*) (27) 

The construction of $(— fe) is essential, but it involves a fair amount of technicality. The 
derivation will be given in appendix B. There, we show that $(— k) is given by 

$(-&) = / Sp U{-p-k,p) 

-^s ■ s r + 7 c r 

5ip(—p) 5?p(p+k) 5ip(—p) 5ip(p + k) 
t K(p/A) 



+ e 



S 



+ 



8^(p)K((p-k)/A) 
K(p/A) 



ip(p — k) 



K((p + k)/A) 



ip(—p — k) 



5 



5t/j(-p) 



S 



(2f 



where the matrix U is defined by 
U {—p — 



= e 



A" ( (p + *)/A) - L^+m^/A) 



(29) 



(30) 



$ is a composite operator, satisfying 

-A-^-^(-k) =X>.$(-Jfe) 

Hence, if (1271) is valid asymptotically for large A 3> /i, then it is valid for any A. As 
explained in sect. 1, this is the most important feature of the Ward identity in ERG. A 
Ward identity similar to (l27j) has been obtained for the scalar QED in equation (7) of 
[TO] , using the background field formalism. 



6. Fine tuning 



We now show how to satisfy the Ward identity (T2~T|) for A ^> fi by fixing the four 
constants (fl3"|) . 
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The asymptotic behavior of k^J^ is easily obtained from the asymptotic behavior 
of the Wilson action, which is given by (I10a|ll0c|ll0a| : 

Sk^J^-k) 



8A v {-k) 



8 



-- ^|a 2 (lnA//i)A 2 

+ 6 2 (lnA//i)m 2 + (c 2 + d 2 ) (lnA//i)fc 2 ) + 



'k^J^y k)- 



8ijj{-p-k)'^"^ ' v <^(p) 
S^M-k) 



(31a) 
(316) 

(31c) 



= a 4 (lnA//i) (kadfjy + kp8 ai + k^ap) + 

The composite operator $(— k) is a dimension 4 scalar operator of momentum — k, 
odd under charge conjugation, and vanishing at k = 0. Hence, its asymptotic behavior 
has the same form as above: 



54,(-A;) _ 

= /c^l^ (In A//j,)A 2 + b 2 (\nA/jj)m 2 + d 2 (\nA/fi)k 2 \ + 
t 



8 



-$(-&)- 



5ip(—p — k) $4>(p) 
5 3 <5>{-k) 



a 3 (lnA / V)^ + 



(32 a) 
{32b) 



(32 c) 



8A a (k 1 )8A /3 (k 2 )8A r (k 3 ) 

= a 4 (lnA//x) (k a 5[3.y + kp8 ai + fc 7 5 ay g) + 

The values of the four dimensionless functions 62, ^2, ^3, 04 at In A/// = determine 
<&(—&) uniquely. Therefore, the operator equation (127)1 is equivalent to the following 
four equations: 

6 2 (0) = 6 2 (0) 
d 2 (0) = J 2 (0) 
Os(0) = o 3 (0) 
04(0) = 04(0) 

where we have used the convention c 2 (0) = 0. 

Now, $ is defined by (I28l 129]) in terms of the Wilson action S(A). Therefore, we 
obtain 



(33) 



8A v {—k) 
t 

8ip(—p — k) 



+ 



/Sp Ui-q-k^jrTZ 



8S 



8 



$(-*;) 



5 



<ty(p) 



e(l-a / (lnA//i))^ 



(34a) 



9 5-0(-p - k) 



Sp < t/(-g- fc,g) 



-S 



8 



Sijj(—q) 8ijj(q + k) f 8ip(p) 



(346) 
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5 3 <5>(-k) 
SAaik^SA^SA^fa) 

= / Sp C/(-g - k,q)- 

J a ' 



5 



5 3 5 



5 



'^(-^^(^^(^^(^^(ff + fe) ( ' 31 ' :) 
Expanding the right-hand sides in powers of the electron mass and external momenta, 
we can obtain b 2 , d 2 , 0,3, 0,4. On the right-hand sides one extra loop is given explicitly 
by the integral over q, and it is sufficient to know 5(A) only up to (/ — l)-loop level to 
determine the four constants 62(0), ^2(0), a 3 (0), 04(0) at /-loop level. (Note a/(0) = 
by convention.) In other words, with 6 2 (0), <i 2 (0), a 3 (0), a 4 (0) determined up to (/ — 1)- 
loop level, we can determine the /-loop values of 6 2 (0), g? 2 (0), a 3 (0), a 4 (0) by demanding 
( 1331) . We can thus construct QED using the Ward identities. This recursive nature of 
the Ward identities has been observed by Becchi [3] and Morris and D'Attanasio [T3] 
among others. We have taken advantage of it using our particular parameterization of 
the theory. 



7. 1-loop calculations 

Using the results of the previous section, let us do 1-loop calculations. We start from 
the tree-level values: 



6 2 (0) = d 2 (0) = a 4 (0) = 0, a 3 (0) 
At 1-loop, (134 ap gives 
<5$ 



8A v {—k) 

Changing the integration variable as 
q — ► Aq 



/ Sp l v U{~q - k,q) 

Jq 



we obtain 



5$ 



8A v {—k) 



e 2 A 2 / Sp lv 

Jq 



K(q) 

A(g)' 
(q 2 ) 2 

where we have expanded the integral in powers of m/A and k^/A. Using 



A J ii + im/A i + tjA + im/A 
K -2A 2 jf lA(g) (1 - K{q)) + [m 2 + jj jf 



A(q)K(q) n 



2\2 



n 



d 4 q -2 d 
(2vr) 4 ^^2 



71+1 



(47r) 2 n 



(independent of the choice of K), we obtain 

2e 2 ... 2e 2 



(35) 



(36) 



(37) 



(3f 



(39) 
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Next we consider (134 fej) at 1-loop: 

- e(l - a/(lnA//i))# 



o 



5ip(—p — k) fi4>(p) 
= ~e 2 / Jt*U(-q - k,q)j„ 



<4 



l-K{{q-p)/A) ( ( g _ p y g _ p ) t 



(q-p) 2 \ 1 " (q-p) 2 J 

~e"i I j^y 2 (1 " K(q)) 2 + ^ (1 - A(g)| (41) 

where we have used the change of variables ( |371) . and expanded the integral in powers 
of k^/A, Pn/A. Hence, 

3-n 



03(0) - e = -e 3 U J^-±- 2 K{q){\ - K(q)) 2 + 



4(4tt) 



(42) 



Finally, we consider (134 eft : 
5 3 ${-k) 



5A a (/c 1 )5A /3 (A; 2 )5A 7 (A;3) 
= e 3 /s P U(-q-k,q) 

J a 



f l-AT((g-fc 3 )/A) l-K((q-k 2 -k 3 )/A) 
\ 7 q-^^ + im k — t-z -iz + im 



Hence, we obtain 



+ (5 permutations) j 

2e 4 (* a <^ + V 7Q + fe^) / (1 - K{q/A)) 2 (43) 

■.«•> - ^ < 44 > 

8. BRST invariance 

To elucidate the gauge structure of QED further, we introduce ghost and antighost fields 
c, c. These fields are non-interacting, and we define the total Wilson action by 

5(A),S(A)-/ t C-(- fc )c W ^ (45) 

We now define the BRST transformation as follows: 

S.A^k) = k^c{k) (46a) 

5 e c(k) = (466) 

S e c(-k) = - ~W(-A;)e (46c) 
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S M - p)= -j k<k)e K ^lf )/K) ^- p - k) («.) 

where e is an arbitrary Grassmann variable. The transformation of ip is highly non- 
linear due to the second term. The Ward identity (|27|) is equivalent to the following 
BRST invariance: 

S,S + I £ c (*) I Sp U( -p - t, p) ^S^^ = (47, 

The first term is the change of the Wilson action S under the transformation. The 
second term is the change of the fermionic integration measure. Note that the non- 
vanishing 5 e S implies the non-invariance of the action under the BRST transformation. 
For example, the non-invariance due to the photon mass term is cancelled precisely by 
the jacobian. 

Due to the electron number conservation, the transformation of ip, ip under BRST is 
not unique. Alternatively, we can define the BRST transformation which is symmetric 
between ip and -0: 

8Mp) =eK(p/A) I ec(k)^(p-k) (48a) 

Jk 

5 e $(-p) = - eK(p/A) f ec(k)^(-p - k) (486) 

Jk 

where 

*<"> ^< p > + 1 ^ A> *X) Si °' (49a) 

•HO = + 3^*1^*1 (496) 

are composite operators. The transformations of A^, c, c are the same as before. The 
BRST invariance is now given as 



5 e S — / Sp 

Jk 



5 5 



(50) 



The first form of BRST invariance has been extended further using the antifield 
formalism in the joint work with Igarashi and Itoh.f 



9. Conclusions 



For some time the compatibility of a smooth cutoff with gauge invariance has been 
established by the previous works ([3],[5],[6],[9]-[T2] and references therein) on ERG. 
Nevertheless, the contrary viewpoint is still prevalent among the practitioners of 
quantum field theory. A main task of this paper was to introduce a concrete application 
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of ERG to convince the reader of the compatibility. In particular we have shown how 
to construct QED using ERG. The Wilson action with a finite momentum cutoff has 
exactly the same gauge invariance as the continuum limit. Nothing is compromised. The 
crucial observation for our construction is that with a Wilson action we can compute 
the correlation functions in the continuum limit. (See (EHH]) and appendix A.) 

As we have mentioned in sect. [U it is some time since the ground work was laid 
down for incorporation of gauge symmetries in the ERG formalism. As far as formalism 
goes, Becchi's results in [3] are seminal, if not final. A merit of the present paper 
is in the use of a particular parameterization of the theory, based on the asymptotic 
behavior of the Wilson action at a very large cutoff. The Ward identities of the Wilson 
action simplify for a large cutoff, and we were able to take advantage of it using our 
parameterization scheme. In a future publication we wish to apply the same scheme to 
YM theories. 
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Appendix A. The cutoff dependence of the correlation functions 

The purpose of this appendix is to explain and derive (J71l8f9~l) and (I18II20I) which play 
important roles in the main text. For simplicity, we derive the analogous results for the 
simpler theory of a real scalar field, which is defined by the Wilson action: 



where 



S(A) = S fae (A) + S„, t (A) (A.r 



I r 

Sint(A) = £ / V 2n (A;p 1 ,---,p 2n )0(p 1 )---0(p 2n ) (A.3) 

^[ [2n)\ J Pl ,-, P 2n 

The interaction part satisfies the following Polchinski differential equation: 

_ aSint(A) = f A(p/A) 1 f 5S int 6S int 6 2 S int \ 

dA J P p 2 + m 2 2\ 5<j){p) 5<t>{-p) 5cj)(p)S(j)(-p) J 1 ' } 

The correlation functions are computed perturbatively using the propagator 

K(p/A) 

p2 _j_ m 2 



(A.5) 
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and the interaction vertices V2 n - We denote the connected part of the n-point correlation 
function as 

(0(pi) ■ • ■ 0(Pn)>5(A) • (2vr) 4 5 (4) (px + ■ ■ ■ + Vn) 

= i/[#]0( Pl )---0(Pn)e 5 ( A ) (A.6) 

where 

Z = J[d(j)]e sw (A.7) 
Similarly, the correlation functions of a composite operator 0(p) is defined by 
(0(p)<f>(pi) ■ ■ ■ 0(Pn)> s( A) • (2vr) 4 5 (4) (p + Pl + ■ ■ ■ + Pn ) 

' r [d<j>}0(p)(f>{p 1 )--- ( f>{p n )e s ^ (A.8) 



Z 

As a preparation, we first prove 

(" A Jv + g Y&M) ■ ■ ■ *<P»))s W = (A.9) 

for any composite operator that satisfies the ERG differential equation 

-A^- = V-0(p) (A.10) 

where 

r A(g/A) / ^ 5 1 g \ 

Ag 2 + m 2 ^0(g)^(-g) 2 50(g)50(-g)y l ' ; 

(This result was first shown in [3].) (Proof) Using ( 1A.4I) and ( 1A.10I) . we obtain 

- A lr {\ f Wi Oip)<t>ipi) ■ ■ ■ <P(Pn)e SW 



dA 

— /[rf0]0(p!)---0(p 



A(g/A) 1 
q q 2 + m? Z 



0(pi) • • • (f)(p n ) 



5<P(-q)5<P(q) 
Afe/A) 1 



5 



2 5<K-q)5<Kq) V ^ ^ 50(g) V" V "^(g/A) 
E §|?^y • § / m O(p)0(Pi) " " " 0(P«> 5 (A. 12) 



which gives (1A.9jl . 

From flA.91) . we find that 

n i 

g ^7X) ' <°(P)^) ■ • • #P») W) (A.13) 

is independent of the cutoff A. 

A particularly important example of a composite operator is 

°o> = ^ (A ' 14) 
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Figure Al. The second graph has at least one interaction vertex. 



We can show that 5S int /5(f)(— p) satisfies (IA.10I) by differentiating the ERG differential 
equation ( 1A.4I) with respect to </>(— p). 

Now, we are ready. Let us first consider the two-point function. Using the Feynman 
rules, we immediately find 



{<t>(p)<t>{-p)) S(A) 



K{p/A) K(p/A) I 5S_ 



P 



+ 



hit 



pi _|_ m z 



\8<K-pY 



'-p) 



(A.15) 



5(A) 



The second term denotes the contribution of all the Feynman diagrams with at least 
one interaction vertex. (See Fig. Al.) This is an analogue of (|T8|) . Hence, we obtain 



K{p/Kf 



(0(P)0(-P))S(A) " 

1 1 / 5S int 



K(p/A) 



p2 _|_ m 2 



'-P) 



K{p/A)p* + m? \do{ P ) 
Since 5Si nt /S(j)(—p) is a composite operator, this is independent of A. Thus, 



(<f>(p)<f>{-p)) 



P 



+ 



<0(P)0(-P)) S (A) 



K{p/ A) 



p2 _|_ m 2 



(A.16) 



(A.17) 



is independent of A. This is the analogue of fl7|5|) . 

For the higher-point correlation functions, there is no graph without interaction 
vertices. (See Fig. A2.) Hence, we obtain 

v u Kipx/k) / 5S int 



-(j)(p 2 ) ■ ■ ■ (j)(p n ) 



(A.18) 



p\ + m 2 \5(f)(-pi) , 5(A) 

for n > 1. This is an analogue of fl20|) . Since SSi nt /5(p(—pi) is a composite operator, we 



find 



5S; r 



pf + m 2 \S(j)(—pi) 

= (0Ol)"--0(Pn)) 



<j)(p 2 ) ■ ■■(f>(p n ) 



n 



S(A) L\ K(pi/A) 



(A.19) 




Figure A2. There is at least one interaction vertex for the higher-point correlation 
functions. 
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is independent of A. This is the analogue of 

An alternative derivation of the results given above can be found, for example, in 
sect. 13 of pQ or in sect. 2 of |14| . where more general functional methods with external 
sources are used. 



Appendix B. Construction of $(— k) 

The composite operator $(— k), introduced in section [51 is given as the sum of two 
operators: 

= di-k) + 2 (-k) 



(b.i; 



where 



Oi(-k) = / S V U(-p-k,p) 

J V 

If 



-S-S- 



5 



+ 



s s- 6 



2 (-k) = e ! 

Jp 



5ip(—p) 5ip(p + k) 5ifj(—p) 5tfj(p + k) 
1 K(p/A) 



(B.2) 



- S 



+ 



<ty(p) K{(p-k)/A) 
K(p/ A) 



ijj{p — k) 



-ip(—p — k)- 



5 



-S 



(B.3) 



K(( P + k)/Ay y c '~ y 5tp{- P y 

Neither 0\ nor 2 are composite operators by themselves. (They do not satisfy 
—A-^Oi = V ■ Oi.) We wish to compute the correlation functions for Oi, 2 in order 
to show that $ has the desired correlation functions (1261) . 

The operator 2 {—k) is the change of the Wilson action under the following linear 
change of variables: 

K(p/A) 



H(p) 

5%[){—p) = f - 



'K((p-k)/A) 
K{p/A) 



ip{p — k) 



ip(—p — k) 



(B.4) 
(B.5) 



■K((p + k)/A) 

Since the correlation functions are invariant under linear changes of variables, the 
correlation functions of 2 {—k) are given by 

'p 2 {-k)A^{k 1 ) ■ ■ ■ A m (k M )ip(pi) • ■■^(PN)^(-qi) ■ ■ 



5(A) 



N 



-E 

i=i 



(A^ (k,)--- 5ij( Pi ) ■ ■ -) 5(A) + (A^ (h) ■ ■ ■ 6$(- qi ) ■ 



5(A) 



(B.6) 



This gives 



TlZi K{h/A) nf =1 K( Pj /A)K( qj /A) 

[0 2 (-k)A fll (k 1 ) ■ ■■A IXM (k M )ip(pi) ■ ■■^(PN)^(-qi) ■ ■•V'(-giv) 

N 

e£ [(A^{h) ■ ■ ■ i)( Pi - k) ■ ■ ■) - (A^{h) ■ -*)••• 
i=i 



5(A) 

(B.7) 
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for arbitrary M, N except for M = 0, N = 1. In the exceptional case we obtain 

' '0 2 {-k)i,{q + k)i,{-q) 



K((q + k)/A)K(q/A)\ v ' rvi yrv *'/s(A) 



e 



K(q/A) 2 \ rv ^ /rv ""/5(A) 



ip(q + k)ip(—q — k) 



K((q + k)/A) 2 \ rv ^ /rv * ; /5(A) 
- k,q) 



K((q + k)/A)K(q/A) 

+ e ((^(q^i-q)) - l^{q + k)${-q - k))) (B.8) 

where we used flSJ) in the last step. The first term on the right-hand side is unwanted, 
and it is the role of 0\ to remove it. 

Let us compute the correlation functions for 0\. For an arbitrary insertion of gauge 
and if), ip fields, we obtain 



f ****** L ishi (• ■ ■ eSW ) ^Lk\ u " ( - q ~ *■ q) = (R9) 



This gives 



TV 



i=l 



^(-Pi + fc) / S(A) 

\ fyfai + k) / S(A) _ 



(B.10) 



The case M = 0, N = 1 is exceptional because of the inhomogeneity of the equations of 
motion: 

(u(-q-k,q)-J^S(A).$(-q)\ = -U(-q-k,q) (B.ll) 
\ 5 ^-P) /5(A) 

L(q + k)-S(A) 5 U(-q-k,q)\ = -U{-q-k,q) (B.12) 
\ Si/j(q + k) / S(A) 

Hence, we obtain 

(Oi(-fe)V(g + *M-?)) 5(A) = "^(-9 " fc > (B.13) 
Thus, with ( IB. 81) we obtain 

[<$>(-k)ip(q+ k)i){-q)\ 



1 '(0 1 + 2 )(-k)ij(q + k)ij(-q) 



K((q + k)/A)K(q/A)\ K 1 z/ v /rv * /rv *'/s(A) 

e ((v(g)^(-g)) - (V>(<? + - fc))) (B-14) 
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For the cases other than M = 0, N = 1, the equations of motion are homogeneous: 

{■■ u{ - p " p -- k) wirT) sw --)s ( ,r (R15) 

{■■ S(A) mTT) U{ - qt - k - qi) --) sw ^° (R16) 

(Note that only the connected parts are considered.) Therefore, we obtain 



OiHOiV (h) ■ ■ • V»(pi) ■ • • ■ ■ ") 5(A) = (B.17) 

Hence, from ( IB. 71) we obtain the desired result 

(H-k)A m (ki)--^(p 1 )---$(-q 1 )---) 

N 

= e [(An ' " -ipiPi — k)-- ■) — (A,, (h) ■ ■ ■ j,(- qi - k) ■ ■ ■)] (B.18) 

i=l 

References 



[l] 
[2] 
[3 

[4] 
[5] 
[G] 
[7] 
[8 
[9 
[10 

[11 
[12 

[13 
[14 



Wilson K and Kogut J 1974 Phys. Rept. C12 75 
Polchinski J 1983 Nucl. Phys. B231 269 

Becchi C 1991 On the construction of renormalized gauge theories using renormalization group 

techniques (Parma lectures) Preprint hep-th/9607188 
Warr B 1988 Ann. Phys. 183 1; Warr B 1988 Ann. Phys. 183 59 
Ellwanger U 1994 Phys. Lett. B335 364 

Bonini M, D'Attanasio M and Marchesini G 1994 Nucl. Phys. B418 81 

Sonoda H 2003 Phys. Rev. D67 065011 

Sonoda H 2007 J. Phys. A40 5733 

Reuter M and Wetterich C 1994 Nucl. Phys. B427 291 

Freire F and Wetterich C 1996 Phys. Lett. B380 337 

Arnone S, Morris T and Rosten O 2005 JEEP 0510 115 

Morris T and Rosten O 2006 Phys. A39 11657 

Morris T and D'Attanasio M 1996 Phys. Lett. B378 213 

Igarashi Y, Itoh K and Sonoda H Quantum master equation for QED in ERG 
Preprint arXiv : 0704 . 2349 



